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ABSTRACT 

The  first-stage,  powered- flight  trajectory  of  a  large  rocket  pow- 
ered vehicle  is  studied  by  varying  the  initial  acceleration,  the  vertical 
flight  time,  and  the  initial  tilt  angle.   Trajectories  were  computed  on  an 
IBM  65O  digital  computer.   Specific  areas  of  interest  with  respect  to  high 
initial  acceleration  rockets  are  the  feasibility  of  using  the  "gravity  turn" 
maneuver  to  obtain  low  burnout  flight  path  angles,  and  the  determination  of 
maximum  energy  trajectories  for  various  values  of  initial  acceleration. 

Results  indicate  that  a  relatively  low  initial  tilt  angle  followed 
by  a  "gravity  turn"  maneuver  is  not  adequate  to  achieve  low  burnout  flight 
path  angles  for  high  initial  acceleration  vehicles.   For  values  of  initial 
acceleration  of  about  2.5  to  3'0  a  large  percentage  of  burning  time  is  spent 
in  the  programmed  tilting  phase,  which  results  in  lift  load  factors  of  the 
order  of  .8  to  1.2. 

Maximum  energy  trajectories  occur  at  specific  values  of  burnout 
flight  path  angle  for  the  initial  accelerations  considered.   These  burnout 
angles  start  at  about  fifty- five  degrees  for  an  initial  acceleration  of  3.0 
and  decrease  to  approximately  zero  degrees  for  an  initial  acceleration  of  1.5 • 

Burnout  conditions  of  velocity,  altitude,  energy,  and  flight  path 
angle  are  plotted  for  the  trajectories  computed.   The  trajectories  most  closely 
approximating  the  maximum  energy  cases  are  included  in  tabular  form. 

Thesis  Supervisor:   Paul  E.  Sandorff 

Title:   Professor  of  Aeronautics  and 
Astronautics 
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OBJECT 


The  object  of  this  thesis  is  to  study  the  early  powered-flight 
trajectory  of  a  large  rocket  powered  vehicle.   The  effects  on  the  first- 
stage  trajectory  of  varying  vertical  flight  time,  initial  tilt  angle, 
and  initial  acceleration,  are  of  primary  interest,  especially  as  they 
affect  the  maximum  burnout  energy  conditions. 

Of  interest  also  is  the  feasibility  of  using  a  relatively  small 
initial  tilt  angle  followed  by  a  "gravity  turn"  to  reach  practical  burnout 
conditions  of  velocity,  altitude,  and  flight  path  angles  for  vehicles  with 
high  initial  acceleration. 


CHAPTER  1 


INTRODUCTION 


The  study  reported  herein  is  concerned  primarily  with  the  initial 
portions  of  the  powered-flight  trajectory  of  a  large,  single  stage,  rocket 
powered  vehicle.   Conceptually,  this  vehicle  could  he  the  "booster  stage  of 
an  ICBM  or  a  satellite  launcher. 

Usually  there  are  three  phases  to  the  initial  flight  trajectory 
of  a  large  ballistic  missile  or  satellite  launching  vehicle.   These  phases 
include  a  vertical  flight  phase,  a  tilt  phase,  and  a  gravity  turn  phase. 
The  gravity  turn  phase  is  customarily  followed  by  a  period  of  "constant- 
attitude  thrust",  during  which  the  major  portion  of  the  flight  velocity  is 
achieved;  this  latter  regime  is  not  considered  in  this  study. 

A  vertical  launch  for  a  large,  rocket- powered  vehicle  of  current 
design  is  necessary  due  to  the  inability  of  the  vehicle  structurally  to  with- 
stand the  transverse  loads  which  would  be  present  during  an  inclined  launch. 
Vertical  or  near  vertical  flight  is  also  necessary  in  order  to  achieve 
altitude.   Usually  the  vertical  flight  path  is  followed  for  a  short  time,  but 
the  time  of  vertical  flight  must  be  carefully  selected  in  order  to  achieve  a 
trajectory  which  minimizes  propellant  expenditure. 

Upon  completion  of  the  vertical  flight  phase,  a  tilting  phase  is 
commenced.   Tilting  is  normally  accomplished  by  deflecting  the  thrust  vector 
of  the  vehicle  to  produce  a  tilting  moment  according  to  some  selected  program; 


this  changes  the  attitude  of  the  vehicle,  and  subsequent  thrusting  changes 
the  velocity  vector.   This  maneuver  is  non-optimum  and  is  best  completed 
quickly;  however,  the  tilt  rate  must  not  be  so  rapid  as  to  exceed  practical 
limitations  of  the  vehicle  control  and  structure.  The  tilting  phase  is  com- 
pleted when  the  vehicle  body  axis  and  the  thrust  vector  are  both  aligned 
with  the  vehicle  velocity  vector. 

The  third  phase  of  the  conventional  trajectory  concerns  the  flight 
regime  where  this  alignment  exists  and  a  relatively  slow  turning  path  follows, 
brought  about  by  the  component  of  gravity  transverse  to  the  flight  path. 
This  phase  hopefully  terminates  at  an  altitude  above  the  sensible  atmosphere, 
and  with  an  attitude  that  matches  the  subsequent  constant- attitude  thrust  re- 
gime in  such  a  manner  that  the  best  overall  trajectory  performance  is  obtained, 

The  important  problem  of  proceeding  from  the  earth's  surface, 
through  the  three  phases  of  the  trajectory  outlined,  to  arrive  at  a  desirable 
altitude,  velocity  and  attitude,  is  complicated  because  of  the  external  forces 
acting  on  the  vehicle.  The  major  forces  affecting  the  vehicle  during  these 
phases  are  thrust,  the  earth's  gravitational  force,  and  the  aerodynamic  forces 
of  lift  and  drag,  which  act  in  a  direction  perpendicular  and  parallel  to  the 
instantaneous  direction  of  flight,  respectively.  Gravitational  and  drag 
forces  acting  on  the  vehicle  result  in  velocity  losses  during  the  flight  and 
thus  detract  from  the  efficiency  of  the  launch.   The  lift  force  may  in  certain 
cases  be  beneficial  in  that  it  may  aid  in  turning  the  vehicle. 

For  a  specific  vehicle,  the  important  trajectory  design  parameters 
are  velocity,  altitude,  and  flight  path  angle  at  burnout.   It  is  only  possible 
however,  to  compute  trajectory  characteristics  by  numerical  integration  of  the 
equations  of  motion  from  specified  initial  conditions.   In  this  paper  numerous 


trajectories  are  developed,  using  vehicle  characteristics  which  are  approx- 
imately representative  of  large  chemical  rockets  of  contemporary  design,  and 
varying  the  time  of  vertical  flight  and  the  maximum  tilt  angle.   The  effects 
of  variations  in  two  important  vehicle  design  parameters  are  also  included: 
namely,  the  initial  thrust-to-weight  ratio,  n^,  and  the  mass  ratio;  the  value 
of  n^  is  introduced  as  an  additional  initial  variable,  while  with  the  assump- 
tion of  constant  mass  flow  every  point  in  each  computed  trajectory  corresponds 
to  burnout  for  some  specific  mass  ratio.   The  burnout  conditions  are  then 
examined  as  functions  of  the  initial  variables  by  using  burnout  angle  as  a 
governing  parameter  and  cross-plotting.   The  nature  of  trajectory  optimization 
to  maximize  burnout  velocity  or  burnout  energy  is  of  particular  concern  to 
this  study. 


CHAPTER  2 


VEHICLE  DESCRIPTION  AND  AERODYNAMICS 


This  study  is  intended  to  derive  conclusions  applicable  to  rocket 
vehicles  similar  to  contemporary  long  range  "ballistic  missiles,  satellite 
launchers,  and  space  vehicle  boosters.   Development  of  trajectory  data  re- 
quires the  use  of  certain  vehicle  design  parameters  which  identify  aerody- 
namic and  engine  performance.   To  simplify  preliminary  work  the  "high-drag" 
configuration  missile  of  Ref .  1  is  selected  as  a  model.   It  is  believed 
that  this  design  has  aerodynamic  characteristics  representative  of  the  class 
of  vehicles  described  above.   Rocket  engine  specific  impulse  is  taken  to  be 
300  lb-sec/lb,  and,  again  for  simplification,  this  value  is  considered  con- 
stant throughout  the  flight  regime.   A  value  of  10  is  selected  for  the 
ratio  of  initial  weight  to  burnout  weight,  defined  as  the  mass  ratio.   This 
makes  the  propellant  factor,  the  ratio  of  initial  fuel  weight  to  initial 
vehicle  weight,  equal  to  .9.   Figure  1  shows  the  physical  dimensions  of  the 
selected  vehicle. 

Since  the  vehicle  of  this  study  is  similar  to  the  high  drag  missile 
of  Ref.  1,  the  aerodynamic  coefficients  utilized  are  extracted  from  this 
source.   A  detailed  explanation  of  the  methods  and  procedures  used  to  arrive 
at  these  values  are  set  forth  in  Appendix  G  of  that  report. 

The  missile,  being  axi ally- symmetric,  has  only  a  drag  force  imposed 
on  it  during  the  vertical  flight  phase  and  the  gravity  turn  phase,  since 
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Fig.  1   Missile  configuration.   (Reference  l) 


during  these  phases  the  angle-of- attack  is  zero.   The  zero-lift  drag  force 
is  made  up  of  three  parts:   base  drag,  skin  friction  drag,  and  form  drag. 
The  zero-lift  drag  coefficient,  based  upon  both  theoretical  and  empirical 
data,  and  representing  the  sum  of  these  forces,  is  plotted  versus  Mach  num- 
ber in  Fig.  2. 

In  order  to  simplify  computer  programming,  the  curve  of  C-p,  versus 
Mach  number  is  divided  into  five  segments.   Each  segment  of  the  curve  is 

then  represented  by  a  straight  line  function.   The  breakdown  of  the  Cp) 

o 

curve  and  the  approximating  straight  line  functions  are  shewn  in  Table  I. 
This  straight  line  approximation  of  the  curve  representing  Cn  ,  while  being 
an  approximation,  is  considered  to  be  sufficiently  accurate  for  the  problem 
at  hand. 

During  the  tilting  phase  of  the  trajectory,  the  missile  is  sub- 
jected to  lift  forces,  as  well  as  drag  forces,  since  the  missile  has  an  angle 
of  attack  during  transition  from  the  vertical  flight  phase  to  the  zero-lift 
phase.   Reference  1  outlines  a  cross-flow  method  of  predicting  lift  and  drag 
on  bodies  of  revolution  at  an  angle  of  attack.   In  this  method  the  flow  over 
the  missile  is  separated  into  two  components:   one  along  the  axial  direction 
of  the  body,  and  one  component  normal  to  the  axis.   The  axial  flow  exerts  a 
force  on  the  body  in  the  axial  direction  while  the  cross  flow  exerts  a  force 
in  the  normal  direction.   Reference  1  derives  equations  for  C-r  and  C^  using 
this  theory  of  cross  and  axial  flow. 
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TABLE  I 

STRAIGHT  LINE  APPROXIMATIONS  OF  THE  ZERO- LIFT  DRAG  COEFFICIENT  CURVE 
FOR  VARIOUS  VALUES  OF  MACH  NUMBER 


M 

% 

o 

0.72 

.130 

1.25 

.85OM  -    A82 

1.90 

.983     -    -323M 

l.k 

.522     -    .080M 

7.3 

.328      -    .023M 

.155 

TABLE  II 

STRAIGHT  LINE  APPROXIMATIONS  OF  THE  CROSS-FLOW  DRAG  COEFFICIENT  CURVE 
FOR  VARIOUS  VALUES  OF  MACH  NUMBER 


M 

\ 

0 

.80 

•  75 

3.23Mc  -   .2*9 

2.36  -   .573Mc 

10 


The  terra  C-p,  is  a  drag  coefficient  due  to  the  cross  flow  component.   A  plot 
uc 

of  C-r,  versus  Mach  number,  assumed  to  apply  for  this  study,  is  shown  in 
Fig.  3-   This  plot  is  a  series  of  straight  line  approximations  of  the  cross- 
flow  drag  characteristics  derived  in  Ref .  1  for  the  missile  configuration  of 
Fig.  1.   These  straight  line  approximations  are  described  by  functions  as 
set  forth  in  Table  II.   The  straight  line  approximations  are  considered 
sufficiently  accurate  since,  as  it  can  be  seen  from  the  above  equation,  the 
effect  of  CD  is  minor  at  small  angles  of  attack. 
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CHAPTER  3 


TRAJECTORY  ANALYSIS 


The  portion  of  the  powered  flight  trajectory  of  interest  in  this 
study  is  considered  in  three  phases,  as  discussed  in  Chapter  1.  The  first 
phase  or  vertical  flight  regime  is  followed  by  the  tilt  phase  during  which 
the  vehicle  is  tilted  from  the  vertical  at  a  rate  of  two  degrees  per  second. 

In  Ref .  1  the  tilt  phase  is  approximated  by  impulsive  tilting  to 
5.5  degrees  from  the  vertical  during  a  one  second  time  interval  at  the  end 
of  vertical  flight  time,  followed  by  a  gravity  turn  which  continues  until 
the  desired  conditions  of  attitude,  altitude,  and  velocity  are  reached.   The 
vehicle  is  assumed  to  be  in  the  gravity  turn  as  soon  as  the  impulsive  tilt- 
ing is  accomplished.   The  impulsive  tilting  during  a  one  second  time  interval 
is  justified  by  determining  that  the  required  vehicle  response  time  is  less 
than  one  second  for  the  5-5  degree  tilt  angle.   This  computation  is  made  on 
the  basis  of  the  time  required  to  tilt  through  the  specified  angle  with  the 
maximum  tilting  moment  available  acting  on  the  moment  of  inertia  of  the 
vehicle.   For  the  present  study,  which  is  concerned  with  higher  values  of  n. 
and  consequent  higher  dynamic  pressures  during  tilting,  a  tilt  rate  of  2 
degrees  per  second  is  selected  as  a  reasonable  maximum  value.   This  is  per- 
haps lower  than  necessary  for  tilting  at  sea  level,  but  to  have  a  basis  for 
comparison,  this  rate  is  used  for  all  trajectories  computed.   In  the  computer 
this  tilt  rate  is  approximated  by  increasing  the  tilt  angle,  U,  two  degrees 
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per  second  until  the  tilt  angle  reaches  the  specified  maximum  programmed 
tilt  angle,  U  .   This  value  of  tilt  angle  is  then  held  constant  until  the 
angle-of -attack  of  the  missile  becomes  zero.   At  this  time  the  tilt  phase 
ends  and  the  zero  angle-of -attack  or  gravity  turn  phase  begins.   During 
the  tilt  phase  thrust  is  considered  to  act  parallel  to  the  vehicle  axis. 
The  component  of  thrust  required  for  tilt  is  considered  a  negligible  loss 
compared  to  the  total  thrust  vector. 

In  the  zero  angle -of- attack  phase  thrust  acts  in  the  direction  of 
the  instantaneous  velocity  vector,  which  is  also  parallel  to  the  missile 
axis.   Turning  is  accomplished  by  the  action  of  the  earth's  gravitational 
field.  This  part  of  the  trajectory  would  logically  be  followed  by  a  con- 
stant attitude  or  a  "linear  with  time"  thrust  program,  depending  on  the 
mission  of  the  vehicle.   In  this  study  the  gravity  turn  is  continued  until 
ninety  percent  of  the  missile  mass  is  consumed.   Since  this  paper  deals  only 
with  single  stage  characteristics,  staging  is  not  considered  and  all  results 
pertain  to  first- stage  values. 


Ik 


CHAPTER  k 


EQUATIONS  OF  MOTION 


The  equations  of  motion  are  developed  using  an  inertial  X,  Y 
coordinate  frame.   This  assumes  a  "flat",  non-rotating  earth,  which  is  a 
good  approximation  during  the  early  powered-flight  phase  of  the  type  of 
rocket  vehicle  considered.   The  gravitational  acceleration  due  to  the  earth 
is  assumed  to  be  constant  during  the  portion  of  the  trajectory  of  interest 
in  this  paper.   This  also  is  a  reasonable  assumption  when  the  altitude 
reached  is  small  compared  with  the  radius  of  the  earth,  as  it  is  in  this 
study. 

Rocket  engine  characteristics  are  simplified  by  assuming  constant 
thrust  and  constant  mass  flow  rate.   Both  of  these  quantities  usually  vary 
with  atmospheric  pressure,  thrust  increasing  and  mass  flow  rate  decreasing 
as  altitude  is  increased.   This  means  that  the  specific  impulse  actually  in- 
creases with  altitude  and  that  the  initial  thrust-to-weight  ratio  is  based 
on  the  lower  level  of  thrust  found  at  sea  level.   The  simplifications  made 
in  this  study  specify  a  constant  specific  impulse  of  300  seconds,  which  may 
be  thought  of  as  representing  an  average  value.   The  initial  thrust-to-weight 
ratio  in  this  study  is  therefore  somewhat  larger  than  it  would  be  for  an 
actual  vehicle  of  comparable  performance. 
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Considering  the  vehicle  as  a  point  mass  the  equations   of  motion  are 
D2Y/DT2  „    (F/m)cos  u  _  gave   _(D/m)sin^    -(L/rrOcostf  (l) 

D2X/DT2  *    (F/m) sin  U  -(D/m)costf    +   (L/ra)sintf  (2) 


wherein 


F/m  =  thrust  per  unit  mass 

Save  "  gravi National  acceleration  due  to  the  earth 

D/ra  =   drag  per  unit  mass 

L/m  =  lift  per  unit  mass 

The  lift  terms  are  considered  positive  in  sign  for  the  negative  angle-of- 
attack  condition  which  occurs  during  the  tilt  phase  of  the  trajectory.   Figure 
h   shows  the  vector  relationships  involved. 

Thrust,  lift,  and  drag  forces  per  unit  mass  are  computed  using  the 
nomenclature  of  Ref .  2,  in  which  Tu  is  defined  as  a  fictitious  time  when  the 
total  mass  of  the  vehicle  would  be  consumed. 

From  the  conventional  relationships  between  rocket  parameters 

Is  =  F/w  (h) 

m  =  W./go(l  -  T/Tu)  (5) 

c  ■  IsgQ  (6) 

Applying  (k) ,    (5),  and  (6)  to  the  various  accelerations  due  to  thrust,  lift, 
and  drag  in  (l)  and  (2)  gives 

F/m  =  c/(Tu  -  T)  (7) 
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Fig.  h       Simplified  diagram  of  vector  quantities  associated  with 
the  missile. 
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l  «„2 


D/m  =    (iev%Ag0Tu)/W.(Tu-T) 


i  o„2 


L/ra  =    (lev^LAg^J/WiCVT) 


(8) 
(9) 
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CHAPTER  5 


COMPUTER  PROGRAM 


The  computer  used  in  this  study  is  an  IBM  65O  digital  computer 
located  in  the  computation  center  of  the  Instrumentation  Laboratory  of  the 
Massachusetts  Institute  of  Technology.   Although  not  comparable  in  speed  to 
the  larger  digital  computers  such  as  the  IBM  70^,  it  is  adequate  for  this 
study,  computing  an  average  trajectory  in  about  ten  minutes.   The  computer 
program  is  prepared  using  the  MAC  programming  system  developed  by  the 
Instrumentation  Laboratory  computation  center. 

Time  intervals  for  integration  are  varied  according  to  the  phase 
of  the  trajectory.   During  the  vertical  flight  phase  the  time  interval  is 
set  at  four  seconds  for  vertical  flight  times  of  four  seconds  and  above,  and 
one  second  for  vertical  flight  times  less  than  four  seconds.   The  time  in- 
terval is  reduced  to  one  second  during  the  tilt  phase  to  maintain  comparable 
accuracy  in  computing  the  rapidly  changing  trajectory  quantities.   At  the 
completion  of  the  tilt  phase  the  time  interval  is  increased  again  to  four 
seconds  and  is  held  constant  until  burnout. 

The  initial  conditions  for  the  equations  of  motion  are  set  equal 
to  zero  for  each  run.   Parameters  held  constant  for  all  runs  are:   S/A,  W^/A, 
DU/DT,  I  ,  and  w.   Variable  parameters  for  each  run  are:   T  ,  Tu,  and  Um. 
The  fictitious  burn-up  time,  Tu,  equals  Is/n^.   Since  Is  is  held  constant, 
T  is  directly  proportional  to  n . .   Table  III  lists  the  numerical  values  of 
the  constants  and  parameters  used. 
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TABLE  III 


NUMERICAL  VALUES  OF  CONSTANTS  AND  PARAMETERS 


Symbol 

Value 

Description 

e0 

.0023769  slug/ft3 

Atmospheric  density  at  sea 
level 

So 

32.17^05  ft/sec2 

Gravitational  conversion 
factor 

gave 

32.0  ft/sec2 

Gravitational  acceleration 
acting  on  vehicle  (assumed 
constant) 

S/A 

11.0 

Ratio  of  planform  area  to  cross 
section  area 

W./A 

3010  lb/ft2 

Ratio  of  initial  weight  to 
cross  section  area 

du/dt 

2  deg/sec 

Tilt  rate 

h 

300  sec 

Average  specific  impulse 

n. 

1 

T 
u 
sec 

sec 

3.0 

2.5 
2.0 

1.5 

100 
120 
150 
200 

90 

1C8 

135 
180 
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Tilt  rate  is  held  constant  at  two  degrees  per  second  until  U  is 
reached.   This  is  mechanized  on  the  computer  by  increasing  U  instantaneously 
at  the  beginning  of  each  one  second  time  interval  until  U  equals  U  .   At 
this  point  Um  is  held  constant  until  the  angle  of  attack  becomes  zero.   Lift 
and  drag  are  computed  during  the  tilt  phase  in  the  manner  shown  in  Chapters 
2  and  k.     These  calculations  are  made  at  the  beginning  of  each  time  interval 
and  are  integrated  as  constants  within  the  differential  equation  loop  of  the 
program.   The  error  introduced  by  this  approximation  was  small,  as  a  result 
of  the  selection  of  integration  intervals;  in  general  the  change  in  aerody- 
namic force  from  one  interval  to  the  next  did  not  exceed  three  percent. 
When  the  angle -of- attack  becomes  zero,  U  is  set  equal  to  (90°  -  ~&  ) ,    and 
thereafter  varies  directly  with  tf  .   This  point  marks  the  beginning  of  the 
gravity  turn  phase. 

In  the  gravity  turn  phase  lift  is  set  equal  to  zero  and  drag  is 
calculated  in  the  same  manner  as  above  using  the  zero  angle-of-attack  drag 
coefficient.  The  program  is  terminated  when  T  equals  .9  T ,  which  corres- 
ponds to  the  mass  ratio  of  ten  mentioned  in  Chapter  2.   Values  of  velocity, 
altitude,  range,  flight  path  angle,  and  energy  are  punched  for  each  computer 
time  interval,  and  for  burnout. 
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CHAPTER  6 


RESULTS 


For  the  study  and  results  as  presented  here,  all  vehicle  and 
trajectory  parameters  are  held  constant,  except  the  initial  thrust-to- 
weight  ratio  of  the  vehicle,  the  time  of  vertical  flight,  and  the  maxi- 
mum programmed  tilt  angle.   The  initial  thrust-to-weight  ratios  used  are 
1.5,  2.0,  2.5,  and  3-0.   The  vertical  flight  times  are  varied  from  1 
second  to  2k   seconds,  and  the  maximum  programmed  tilt  angle  is  varied 
from  2  to  90  degrees.   The  trajectory  calculations  are  continued  in  all 
cases  for  a  total  time,  T  =  0.9TU,  i.e.,  a  mass  ratio  of  10. 

Approximately  one  hundred  trajectories  were  computed  for  this 
paper.   Table  IV  lists  values  of  mass  ratio,  velocity,  altitude,  flight 
path  angle,  energy,  scalar  velocity  loss  due  to  drag,  and  scalar  velocity 
loss  due  to  gravity  for  some  of  the  more  useful  trajectories.  The  effects 
of  varying  the  parameters,  ni,  Ty,  and  Um,  are  shown  in  Figs.  5  through  18, 
which  display  burnout  values  of  velocity,  altitude,  energy,  and  flight 
path  angle  plotted  against  maximum  programmed  tilt  angle  for  a  mass  ratio 
of  10.   Separate  plots  are  shown  for  each  vertical  flight  time  used. 

Scalar  velocity  loss  due  to  drag  is  shown  in  Fig.5.19  through  22. 
V-pj,  divided  by  the  ballistic  coefficient,  W./C^  A,  is  plotted  versus  burnout 
flight  path  angle  for  each  n-  and  each  Tv.   A  more  general  presentation  is 
obtained  with  the  ballistic  coefficient,  which  uses  the  value  of  Cn  at  Mach 
2.0. 
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Since   side  loading  is   an  important   consideration  in  large  rocket 

design,   plots  of  the  maximum  lift   load  factor  versus  U     for  each  n.    and  each 
°   '  mi 

T'v  are  included  as  Figs.    23  through  26. 

The   length  of  time  that  the  rocket  is   subjected  to  lift   loads   is 
also  of  interest.      Figure  27   shows  the  time  required  to  tilt  the  missile   so 
that  it  will  attain  a  burnout  flight  path  angle  of  thirty  degrees.      This   is 
plotted  against  n^  to  show  the  large   increase  in  time  required  for  tilting 

as  n.    is  increased. 

i 

An  optimization  study  is  made,  based  on  finding  the  combination  of 
parameters  which  would  give  the  highest  specific  energy  at  burnout  for  each 
n.  investigated.   Energy  is  first  maximized  for  fixed  values  of  ft-u  by  plot- 
ting energy  versus  the  value  of  U  corresponding  to  particular  vertical 
flight  times.   The  maximum  energy  points  are  then  cross-plotted  against  ~8  ^ 
and  the  value  of  U  corresponding  to  the  maximum  energy  point  for  the  fixed 
values  of  tf-u.   As  a  cross-check,  this  procedure  is  reversed  so  that  energy 
is  maximized  for  fixed  values  of  U  ,  and  cross-plotted  in  the  same  manner. 
Also  included  in  these  figures  are  the  values  of  burnout  velocity  and  burn- 
out altitude  which  occur  at  the  maximum  energy  points.   The  optimization 
results  are  shown  in  Figs.  29  through  31. 

Representative  trajectories  for  various  values  of  n.  are  identified 
in  detail  in  Table  V.  These  particular  trajectories  were  chosen  because  they 
were  closest  to  the  maximum  energy  trajectories  for  each  case. 
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TABLE  IV 

TRAJECTORIES  INVESTIGATED  SHOWING  BURNOUT 
VALUES  FOR  VARIOUS  VALUES  OF  MASS  RATIO 


] 

n.  s  1-5 
1 

Ty 

3l 

..\ 

MR 

\ 

\ 

A 

-8 
E^xlO 

VD 

VG 

8 

2 

13* 
lAl 

l6l 
180 

1.07 
3.39 
5.13 
10.0 

»  232 

7,190 

10,684 
16,700 

1,463 
285,739 
423,565 
616,000 

87.8 
53.5 
50.9 
49.0 

.001 

.350 

.707 

1.590 

•59 
422 
422 
422 

420 

4,188 

4,694 

5,078 

8 

6 

14 
l42 
162 
180 

1.08 
3^5 
5.26 
10.0 

252 

8,169 

12,086 
18,130 

1,703 
165,530 
207,113 
250,500 

83.6 

14.8 

ll.l 

8.9 

.001 

.387 

.797 

1.670 

•  78 
749 
779 
798 

483 
3,032 

3,155 
3,272 

8 

12 

17" 
1*1 

161 

1.09 
3.39 
5.13 

31^ 
4,843 
2,970 

2,532 
363,871 
140,962 

78.2 
-11.3 
-20.3 

.001 

.129 

.048 

1.6 
5,083 
11,142 

516 
1,874 
1,688 

16 

4 

26* 
lk2 

162 
180 

1.15 
3.45 
5.26 

10.0 

506 
7,204 

10,731 
16, 4oo 

6,21k 

305,300 
458,176 
650,700 

86.0 
61.4 
59.6 
58.O 

.003 

.358 

.723 

1.553 

5.91 
407 

407 
407 

838 

M39 
4,882 

5,393 

16 

8 

27* 
139 
159 
180 

1.16 
3.28 

4.88 
10.0 

529 

7,136 

10,565 

17,100 

6,715 

248,475 
358,888 
525,000 

82.0 
41.8 
38.4 
36.0 

.oo4 

.334 

.674 

1.628 

6.80 
465 
465 
465 

894 

3,899 
4,220 

M35 

16 

12 

29* 
l4l 

161 
180 

1.17 
3.39 
5.13 
10.0 

577 

7,75^ 

11,484 

17,700 

7,782 
212,551 
290,231 
388,000 

78.O 
27.0 

23.3 
21.0 

.004 

.369 

.753 

1.680 

8.79 
556 
559 
559 

949 
3,490 
3,757 
3,9^1 

16 

18 

31* 

139 
159 
180 

1.18 
3.28 

if. 88 
10.0 

627 

7,727 

11,392 

18,181 

8,900 

147,717 
180,348 
216,304 

72.8 
12.9 

9.0 
5.8 

.005 

.346 

.706 

1.722 

11.3 
821 
878 
937 

952 
2,952 
2,980 
3,082 

*  Completion  of  tilt  phase.  (  °<-  =  0). 
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TABLE 

IV  (Conti: 

aued) 

n±   =  1.5 

T 

V 

U 

m 

\ 

MR 

J^ 

Yb 

4 

v10"8 

i 

VG 

2k 

6 

39': 
139 
159 
180 

1.2k 

3.28 
^. 88 
10.0 

82U 

6,755 
10,011* 
16,300 

ll*,793 
291,072 

^39,719 
676,000 

83.6 
66.0 
63.6 
62.3 

.008 

.321 

.61*3 

1.51*6 

21.5 

1*01* 

1*01* 
1*01* 

1,230 
k,3kl 
M32 
5,!+96 

2k 

12 

1*1* 

l6l 
180 

1.26 
3.39 
5.13 
10.0 

878 

7,310 

10,859 
16,830 

16,392 
270,316 
395,063 
566,000 

78.0 

1*6.5 

^3.5 
1*1.5 

.009 

.35^ 

•  717 

1.597 

26.5 

1*1*3 
1*1*3 
1*1*3 

1,326 

3,7^7 
1*,1*98 
1*,927 

2k 

2k 

^5* 
lUl 
161 
180 

1.29 
3.39 
5-13 
10.0 

988 

7,81*9 

11,630 

17,850 

19,658 
190,631 
252,1*72 
326,500 

66.0 

21.7 
17.8 

15.5 

.Oil 

.369 

•  757 

1.695 

1*2.9 
626 
631* 
636 

1.1*29 
3,025 
3,536 
3,711* 

2k 

32 

1*9* 
ll*l 
161 
180 

1.33 
3-39 
5-13 
10.0 

1,107 

7,996 

11,815 

18,095 

23,086 
11*1,693 
168,1*29 
192, k8k 

58.0 

10.8 

6.9 

k.l 

.011* 

.365 

•  752 

1.699 

66.9 
891 

979 
1,092 

1,576 
2,613 
3,006 
3,013 

*  Completion  of  tilt  phase   (<^<=  0) 
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TABLE  IV  (Continued) 


n±   -   2.0 


Tv 

um 

Tb 

MR 

Zl 

A 

*b 

I^xlO"8 

3l 

VG 

1 

4 

4* 

104 
120 
135 

1.03 
3-26 
5.0 
10.0 

133 

7,785 

11,520 

17,853 

263 
235,616 
351,320 
510,742 

86.0 
51-2 
49.1 
47.5 

.0002 
.379 
.777 

1.758 

.05 

658 

659 
659 

152 

2,957 
3,341 
3,688 

1 

10 

7" 
103 
119 
135 

1.05 
3. OP 

if. 83 

10.  u 

238 

7,861 
11,626 

l8,4o6 

811 

175, 494 
247,918 
350,294 

80.2 
30.4 

27.6 
25.3 

.0005 
.365 

.756 
1.807 

.34 

892 
902 
903 

235 
2,047 
2,652 
2,891 

1 

20 

12-x- 

io4 

120 

135 

1.09 
3.26 

5-0 
10.0 

427 

7,979 
11,871 
18,300 

2,392 
125,765 
162,216 
200,000 

69.I 
16.0 

12.9 
11.0 

.002 

•  359 

•  757 
1.748 

2.0 

1,419 

1,532 

1,623 

4oi 
2,002 
2,117 
2,277 

1 

24 

Ik* 

106 
118 
135 

1.10 

3.40 
4.68 
10.0 

507 

8,216 

11,042 

17,980 

3,250 
114,092 
133,796 
165,825 

66.9 

12.0 

9.6 

6.9 

.002 

•  374 

.653 

I.669 

3.4 
1,720 
1,901 
2,203 

4io 
1,844 
1,937 
2,017 

8 

10 

20* 
104 

120 
135 

1.15 
3.26 
5.0 
10.0 

732 
7,670 

11,357 
17,620 

6,989 
251,704 
380,350 
559,000 

80.0 
60.1 

58.9 
58.0 

.007 

•  375 

•  767 
1.729 

16.5 
616 
616 
616 

602 
3,114 
3,547 
3,964 

8 

20 

25* 
105 
121 
135 

1.20 

3.33 
5.16 
10.0 

945 
8,085 

11,998 
18,040 

10,914 
217,596 
318,203 
443,500 

69.8 
4l .  5 

39.1 
38.6 

.008 

.397 

.822 

1 .  762 

22.3 
738 
739 
739 

793 
2,797 
3,113 
3,421 

8 

30 

30* 
106 
118 
135 

1.25 
3.  to 

4.68 
10.0 

1,164 
8,422 

11,317 
18,350 

15,467 
179,763 
231,661 
332,000 

59.8 
28.1 
26.0 
23.9 

.012 
.413 

.715 

1.786 

54.2 
938 
948 
950 

937 
2,420 
2,615 
2,900 

8 

50 

4o* 
104 
120 
135 

1.36 
3.26 
5.0 
10.0 

1,600 

7,666 

11,280 

17,353 

25,169 
101,365 
121,199 
140,363 

40.9 

10.3 

7.0 

4.5 

.021 

.326 

.675 

1.551 

223 
1,832 

2,275 
2,848 

1,142 
1,902 
1,965 
1,999 

*  Completion  of  tilt  phase.   ( c< =  0) 
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TABLE  IV  (Continued) 


r^  ■  2.0 


T 

um 

Tb 

MR 

Vb 

Yb 

4 

E^xlO"8 

Zl 

VG 

16 

16 

37* 
105 
121 
135 

1.31 
3-33 
5.16 
10.0 

1,1+22 

7,8^9 
11,655 
17,715 

2^,960 
255,783 
386,0^5 
557,000 

7I+.I 
59-*+ 
57.8 

56.5 

.018 

.390 

.803 

1.700 

132 
625 
625 
625 

1,051 
3,lJ+6 
3,570 
3,860 

16 

36 

k6* 
106 
118 
135 

l.kk 

3.  ko 

k.6Q 
10.0 

1,782 

8,327 

11,195 

18,180 

36,270 

196,701+ 
257,1+36 
376,000 

5I+.I 
33.3 
31.1+ 
29.0 

.028 

.1+10 

.710 
1.795 

328 
833 
837 
838 

2,1+05 
2,620 
2,81+8 
3,182 

16 

kQ 

51* 
103 
119 
135 

1.52 
3.08 
1+.83 
10.0 

2,030 

7,807 
11,609 
18, kkk 

1+1,966 
11+7,81+3 
200,871 
273,ll+l 

1+2.6 
22.7 
19.7 
17.3 

.031+ 

.352 

•  739 

1.790 

1+1+7 
1,031 
1,068 
1,081 

1,583 
1,962 
2,503 
2,675 

16 

60 

58* 
106 
118 
135 

1.63 
3A0 
Ik  68 
10.0 

2,k63 

8,253 

11,107 

18,102 

1+2,291 
118,317 
139,002 
173,138 

30.1 

12.1+ 

10.0 

7.1+ 

.01+6 

•  379 

.661 

1.690 

601+ 
1,1+36 
1,587 
1,828 

1,663 
2,091 
2,186 
2,270 

2k 

12 

1+9* 
105 

121 
135 

1.1+8 

3.33 
5.16 
10.0 

1,905 

7,7^5 

11,50)+ 

17,1+20 

1^,785 
271,198 
1+13,801+ 
597,ooo 

78.1 
71.2 
70.1 
69.3 

.033 

.387 

.795 

1.759 

3*+7 
590 
590 
590 

1,533 
3,285 
3,756 
M90 

2k 

30 

56* 
10U 
120 
135 

1.59 
3.26 
5.0 
10.0 

2,281 

7,757 

ll,5H 

17,863 

56,032 
226,562 
336,1+81+ 
1+87,562 

60.0 
k8.2 

1+6.0 

1+1+.3 

.01+1+ 

.371+ 

.771 

1.750 

1+68 
670 
671 
671 

1,726 
2,973 
3,338 
3,666 

2k 

50 

65* 
105 

121 

135 

1.76 

3.33 

5.16 

10.0 

2,905 

8,136 

12,11+0 

18,317 

68,81+6 

176,979 
21+7,291+ 
332,589 

lK).2 
28.1+ 
25.6 
23.6 

.061+ 

.388 

.816 

1.780 

626 
81+2 
853 
85I+ 

1,929 
2f6k2 
2,857 
3,029 

2k 

70 

77* 
105 

121 
135 

2.05 

3.33 

5.16 

10.0 

1+,012 

8,105 

12,099 

18,222 

80,808 
121,801+ 

ll+8,53l+ 
175,681+ 

20. 0 

12.1 

9.1 

6.9 

.106 

.367 
.780 

1.717 

8I+5 
1,190 
1,350 
1,531 

2,073 
2,325 
2,1+01 
2,1+1+7 

*  Completion  of  tilt  phase.  (<=*!.=  0) 
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TABLE  IV  (Continued) 


ni  "  2-> 

T 

%_ 

Jk 

MR 

vb 

\ 

xb 

E^xlO"8 

JSl 

VG 

1 

10 

8* 
84 
96 

108 

1.07 

3-33 

5.0 

10.0 

4io 

8,343 

11,927 

18,292 

3,634 
223,832 
326,512 
476,217 

19.7 
59.8 
58.6 
57-6 

.003 

.420 

.816 

1.826 

3.5 
794 
795 
795 

242 
2,483 
2,798 
3,113 

1 

20 

15" 
83 
95 

108 

1.14 

3.24 
4.78 
10.0 

810 

8,200 

11,711 

18,511 

5,638 
189,882 
272,461 
402,888 

69.9 
45.9 
44.2 

42.7 

.005 
.397 
•  773 

1.843 

9.32 
926 
931 
931 

445 
2,224 
2,4l8 
2,758 

1 

4o 

-St 

27 
83 
95 

108 

1.29 
3.24 
4.78 
10.0 

1,527 

8,125 

11,696 

18,595 

17,249 
136,994 
187,321 
264,277 

50.4 
27.3 
25.1 
23.2 

.017 

•  374 

.744 

1.814 

138 

1399 

1447 
1465 

795 
1,826 

1,917 
2,140 

1 

6o 

39* 
83 
95 

108 

1.48 
3.24 
4.78 

10.0 

2,218 

7,450 

io,547 

16,722 

30,093 

86,735 

105,551 

129,978 

30.5 
12.4 

9.8 
7.5 

.034 

.305 

.590 

1.440 

546 
2419 
3038 
3888 

1,016 
1,474 

1,485 
1,590 

8 

20 

29* 
85 
97 

108 

1.32 
3.42 
5.18 
10.0 

1,584 

8,556 

12,286 

18,200 

21,891 
225,510 
328,567 
465,000 

70.0 
57.7 
56.5 
56.0 

.020 

.438 

.860 

1.670 

167 
834 

835 
835 

929 
2,490 

2,779 
3,165 

8 

30 

34* 
82 
94 

108 

i.4o 
3.16 
4.58 
10.0 

1,846 

7,898 

11,284 

18,437 

28,971 
182,871 
263,301 
402,241 

60.2 
46.5 
44.8 
43.2 

.026 

.371 

.721 

1.829 

310 
94l 
947 
947 

1,084 
2,261 
2,449 
2,816 

8 

50 

44* 

84 

96 

108 

1.58 

3.33 
5.0 

10.0 

2,497 

8,367 

12,074 

18,593 

42,862 
143,370 
195,898 
269,376 

40.5 
27.6 
25.5 
23.7 

.045 
.396 
•  792 

1.815 

623 
1284 
1324 
1338 

1,300 
1,969 

2,122 
2,269 

8 

6o 

50* 
82 
94 

108 

1.72 
3.16 
4.58 
10.0 

3,012 

7,750 

11,150 

18,300 

50,266 
112,553 
146,717 
201,300 

30.1 
19.8 
17.5 
16.0 

.062 

.336 

.669 

1.740 

804 
1554 
1700 
1823 

1,414 

1,806 
1,830 
2,077 

*  Completion  of  tilt  phase.   (  <=<=  0) 
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TABLE  IV  (Continued) 
n±   *  2.5 


T 

Um 

Tb 

MR 

A 

Yb 

k. 

E^xlO"8 

V, 

\ 

16 

20 

1+3* 

83 

95 

108 

1.56 
3.2U 
1+.78 
10.0 

2,1+05 
8,038 

11,1+79 
18,300 

1+8,1+61 
220,1+66 
323,1+66 
1+89,000 

70.1 
61+. 0 
62.8 
62.0 

.01+1+ 

.39^ 

.763 

1.811+ 

510 
780 
781 
781 

1,365 
2,532 
2,800 
3,119 

16 

1+8 

56* 
81+ 
96 

108 

1.88 

3.33 

5.0 

10.0 

3,593 

8,385 

12,071 

18,553 

72,691+ 
169,551+ 

237,219 
333,571 

1+2.3 
35.6 
33.6 
32.1+ 

0.088 
.1+06 
.805 

1.828 

801 

1011 

1021+ 
1026 

1,706 
2,221+ 
2,1+25 
2,621 

16 

60 

67* 
83 
95 

108 

2.26 
3.21+ 
4.78 

10.0 

5,030 

8,123 

11,706 

18,611+ 

92,937 
139,010 

186,315 
258,262 

30.0 
25.7 
23.5 
21.5 

0.156 
.375 
•  7^5 

1.816 

1017 
1163 
1209 
1228 

1,823 
2,061+ 
2,11+5 
2,1+58 

16 

70 

70* 
82 

108 

2.1+0 

3.16 
4.58 

10.0 

5,^67 

7,812 

11,21+0 

18,1+71+ 

91,502 
116,156 
11+6,665 
19^,613 

20.2 
17.8 
15.1+ 
13.2 

•  179 

.3^2 

.679 

1.769 

lll+O 

1323 

i'+5i+ 
1576 

1,81+3 
1,975 
1,986 
2,150 

2k 

10 

55* 
83 
95 

108 

1.85 
3.21+ 
ij.,78 
10.0 

3,530 
8,015 

11, ^23 
18,101 

8l+,  985 
236,080 

3^8,703 
530,669 

80.0 
78.1+ 
77-9 
77.5 

.090 

•  397 

.765 

1.809 

629 
703 
703 
703 

1,781 
2,632 
2,93'+ 
3,396 

2k 

20 

59* 
83 
95 

108 

1-97 
3.21+ 
1+.78 
10.0 

3,97^ 
8,031+ 

11,1+61+ 
18,120 

96,975 
226,01+2 
332,126 
502,000 

70.0 

67.5 
66.6 

65.5 

.110 

.395 

.761+ 

1.811 

680 
732 
733 
733 

1,896 
2,581+ 
2,863 
3,31+7 

2k 

ko 

68* 
81+ 
96 

108 

2.3 
3-33 
5.0 
10.0 

5,17^ 

8,312 

11,91+6 

18,361+ 

121,1+21 
200,512 
287,122 
1+12,507 

1+9.8 
V7.7 
1+6.1 
1+5.0 

.173 

.1+10 

.806 

1.819 

816 

8I+1+ 
81+7 
8I+7 

2,060 
2,1+61+ 
2,927 
2,987 

2k 

6o 

83* 

95 

108 

3.21+ 
4.78 
10.0 

8,051+ 
11,636 
18,525 

160,51+8 
216,651+ 
303,506 

30.0 
28.5 
26.6 

.376 

.71+7 

1.813 

983 
1001 
1005 

2,323 
2,1+23 
2,670 

*  Completion  of  tilt  phase.   (  <=<=  0) 
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TABLE  IV  (Continued) 


n.  =  3-0 
1 


T 

V 

um 

A 

MR 

\_ 

\_ 

4 

E^xlO-8 

VD 

VG 

1 

16 

14* 

70 
78 
90 

1.16 
3-33 
4.53 
10.0 

999 

8,603 

11,370 

18,648 

6,606 
197,650 
266,817 
417,702 

74.2 
61.5 

60.7 
60.0 

0.007 
0.434 
0.732 
1.873 

14.1 
942 
944 
945 

417 
2,078 
2,306 
2,607 

1 

24 

20* 
68 
80 
90 

1.25 
3.12 
5.0 
10.0 

1,448 

8,031 

12,257 

18,700 

13,350 
168,876 
262,685 
378,000 

66.1 
53.0 
51.5 
50.5 

0.015 

0.377 
O.836 
1.867 

84.8 
1,051 
1,061 
l,06l 

622 
1,898 

2,202 
2,439 

1 

4o 

30* 
70 
78 
90 

1.43 
3.33 
4.55 
10.0 

2,l4l 

8,519 

11,335 

18,714 

27,467 
148,091 
195,038 
296,016 

50.1 

37.9 
36.6 
35.0 

0.032 
o.4io 
0.705 
1.846 

418 
1,360 
1,384 
1,393 

901 

1,741 
1,901 
2,093 

1 

52 

37* 
69 
81 
90 

1.59 
3.23 
5.26 
10.0 

2,716 
8,072 

12,515 
18,500 

376,444 

118,215 

171,985 
230,000 

38.1 
28.2 
26.1 
25.O 

0.049 
0.364 
O.838 
1.770 

714 
1,665 
1,776 
1,818 

1,050 
1,583 
1,739 
1,882 

8 

10 

27* 
71 
79 
90 

1.37 
3.44 
4.76 
10.0 

1,895 

8,816 

11,682 

18,503 

25,477 
219,839 
298,606 
454,492 

80.1 
76.O 
75.6 
75.1 

0.026 

0.459 
0.778 
1.858 

279 
895 
896 
896 

866 
2,209 
2,472 
2,801 

8 

20 

32* 
68 
80 
90 

1.47 

3.13 

5.0 

10.0 

2,257 

7,975 

12,163 

18,570 

34,766 
184,230 
290,286 
423,328 

70.1 
64.2 
63.2 
62.4 

0.037 
0.377 
0.833 
1.860 

451 
954 
958 
959 

1,012 
2,051 
2,399 
2,671 

8 

1+0 

42* 
70 
78 
90 

1.73 
3-33 
4.55 
10.0 

3,212 

8,559 
11,366 
18,718 

54,483 
163,808 
217,461 
333,471 

50.1 
44.0 
42.7 
4l.5 

O.069 
0.419 
0.716 
1.859 

722 
1,144 
1,156 
1,160 

1,366 

1,917 
2,098 
2,322 

8 

6o 

53* 
69 
81 
90 

2.13 
3.23 
5.26 
10.0 

4,681 

8,135 

12,599 

18,644 

73,877 
120,154 
171,518 
226,562 

30.2 
26.8 
24.6 
23.2 

0.133 
0.369 
0.849 
1.811 

1,148 
1,473 
1,572 
1,607 

1,471 
1,712 
1,859 
1,949 

Completion  of  tilt  phase.   (=<.=  0) 
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TABLE  IV   (Continued) 


n.    =  3.0 


T 

m 

Tb 

MR 

vb 

Yb 

<, 

^xlO-8 

Zl. 

VG 

16 

16 

45* 
69 
81 
90 

1.82 
3.23 
5.26 
10.0 

3,606 

8,261+ 

12,623 

18,550 

71,902 
201,053 
317,617 
451,000 

74.1 
72.0 

71.3 
71.0 

0.088 
0.406 

0.899 
I.865 

728 
871 
874 
874 

1,446 
2,235 
2,533 
2,776 

16 

32 

52* 
68 
80 
90 

2.08 

3.13 

5.0 

10.0 

4,584 

7,986 

12,201 

18,600 

92,482 
174,404 
271,637 
394,000 

58.2 
56.O 
54.4 
54.0 

0.135 
0.375 
0.832 

1.857 

878 
962 

969 
970 

1,598 
1,532 
2,350 
2,63b 

16 

60 

72* 
80 
90 

3.33 
3.57 
5.0 
10.6 

8,487 

9,112 
12,209 
18,735 

144,099 
152,911 
194,233 
265,485 

30.4 
30.1 

29.3 
27.8 

0.406 
0.464 
0.808 
1.840 

1,196 

1,204 

1,230 

1,243 

1,937 
1,984 
2,081 
2,222 

16 

90 

70 
80 
90 

3.33 

5.0 

10.0 

8,139 
11,824 
18,169 

120,218 
134,099 
143,873 

12.1 
6.2 
2.7 

0.370 
0.742 
1.697 

1,431 
1,574 
1,882 

2,050 
2,122 
2,149 

24 

20 

67* 
71 
79 
90 

3.03 
3.44 

k.16 
10.0 

7,747 

8,871 

11,746 

18,550 

186,250 

217,397 
293,980 
442,800 

70.1 
70.0 
69.4 
69.O 

O.360 
0.463 
0.784 
I.863 

827 
830 
831 
831 

2,126 
2,219 
2,473 
2,819 

2k 

ko 

70 

80* 
90 

3.33 

5.0 

10.0 

8,499 
12,l6l 
18,550 

186,774 
265,425 
379,300 

51.5 
50.0 

49.5 

0.421 
O.825 
1.842 

949 
953 
954 

2,172 
2,4o6 
2,696 

2k 

6o 

70 
80 
90 

3.33 

5.0 

10.0 

8,317 
12,019 
18,528 

3-69,663 
227,544 
327,434 

37.8 
34.0 
32.0 

0.400 

0.795 
1.808 

1,047 
1,057 
1,059 

2,256 

2,444 
2,613 

2k 

90 

70 
80 
90 

3.33 

5.0 

10.0 

7,968 
11,449 
17,876 

162,805 
3-98,428 
230,643 

29.3 
18.1 
10.5 

0.370 
0.719 
1.672 

1,094 
1,119 
1,135 

2,558 
2,952 
3,189 

*  Completion  of  tilt  phase.      (  °C-  0) 
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Fig.  5   Variation  of  V^  with  Um  for  various  values  of  n.  at  a 
T     of  1  and  8  seconds. 
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Fig.  6   Variation  of  Vh  with  U  for  various  values  of  n.  at  a 


Tv  of  16  and  2H-   seconds. 
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Fig.  7   Variation  of  Y,  with  U  for  various  values  of  n.  at  a 
T  of  1  second?      m  1 
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Fig.  8   Variation  of  Y,  with  U  for  various  values  of  n  at  a 


T  of  8  seconds. 
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Fig.  9   Variation  of  Y^  vith  U  for  various  values  of  n.  at  a 

T  of  16  seconds.  1 
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Fig.  10   Variation  of  Y^  with  Ura  for  various  values  of  n.  at  a 
T  of  24  seconds. 
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Fig.  11   Variation  of  K  with  U  for  various  values  of  n  at 
T  of  1  second.  ± 
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Fig.  12   Variation  of  E.  with  U  for  various  values  of  n  at  a 
T  of  o  seconds. 
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Fig.  13   Variation  of  F^  vith  Um  for  various  values  of  n.  at  a  T 
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Fig.  Ik       Variation  of  E,  with  U  for  various  values  of  n.  at  a  T 
of  24-  seconds. 
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Fig.  15   Variation  of  ^f  h  with  U_  for  various  values  of  n,  at 


Tv  of  1  second. 
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Fig.  16   Variation  of   )fb  with  Um  for  various  values  of  n.  at  a 


Tv  of  8  seconds. 
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Fig.  18   Variation  of  0^  with  Um  for  various  values  of  n. 
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Fig.  19   Drag  velocity  loss  as  a  function  of  #*,  for  various 
values  of  n±   at  a  Tv  of  1  second. 
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Fig.  20   Drag  velocity  loss  as  a  function  of  X     for  various 

values  of  n,  at  a  T  of  8  seconds, 
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Fig.  21   Drag  velocity  loss  as  a  function  of  0,  for  various 


values  of  n.  at  a  T  of  l6  seconds, 
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Fig.  22   Drag  velocity  loss  as  a  function  of  #   for  various 

values  of  n.  at  a  T  of  2k   seconds. 
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Fig.  23   Variation  of  maximum  lift  load  factor  with  U  for  various 
values  of  n.  at  a  T  of  8  seconds. 
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Fig.  2k      Variation  of  maximum  lift  load  factor  with  U  for  various 


values  of  ^  at  a  T  of  8  seconds. 
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Fig.  25   Variation  of  maximum  lift  load  factor  with  U  for  various 


values  of  n^.  at  a  T  of  l6  seconds. 
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Fig.  26   Variation  of  maximum  lift  load  factor  with  U  for  various 

values  of  n.  at  a  T  of  2k   seconds.        m 
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Fig.  27   Variation  of  tilting  time  with  n.  for  various  values  of  T 
to  reach  a  fr,  of  30°. 
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under  maximum  energy  conditions  for  n.  of  2.0. 
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TABLE  V 


COMPUTER  RESULTS  FOR  REPRESENTATIVE  TRAJECTORIES 


A. 

n.  =  3-0 

1 

T  =  1  sec 

U  =  20° 
m 

T 

V 

Y 

If 

E 

VD 

sec 

fps 

ft 

dee; 

ft-lb/slUR 

fps 

4 

265.78352 

525.70761 

87.002181 

52234.583 

. 20644683 

8 

5^6.99767 

2129.7715 

80.949374 

2l8l26.6l 

2.2177488 

16 

1152. 0408 

8609. U335 

70.776125 

940599.34 

27.850098 

17* 

1227.7526 

9726 . 6649 

70.160493 

1066634.4 

37.685454 

21 

1535.2739 

11*867.757 

69.635561 

1656889.0 

87.560342 

29 

2087.3356 

28129.005 

66.455353 

3083509.0 

331.30749 

37 

2731.9199 

45316. 107 

63.757048 

5189695.9 

610.93239 

4? 

3626.1614 

67522.151 

61.496942 

8746984.1 

802.44887 

53 

4799.5585 

96429.217 

59-667247 

14620399-0 

925.23353 

6l 

6325.6918 

13^027.09 

58.208353 

24319383.0 

982.43218 

69 

8308.6788 

182846 A5 

57.055556 

40399983.0 

1000 . 4780 

77 

10972.134 

246406.84 

56.151720 

68121766.O 

1005 . 5098 

85 

14886.326 

330622.30 

55.454918 

121438810 

1006.2890 

90 

18668.867 

398806.33 

54.922312 

187094520 

1006.2980 

*  End  of  tilt  phase. 


B. 

n.  ■  2.5 

Tv  =  1  sec 

um  ■  26° 

T 

V 

Y 

1 

E 

VD 

sec 

fps 

ft 

de£ 

ft-lb/slug 

fps 

4 

199.15440 

393-99420 

86.681749 

32507.627 

.11541601 

8 

409.79788 

1593-1459 

80.054237 

i35225.ll 

1.2274019 

12 

632.93575 

3600.1614 

73.392714 

316135.60 

4.6832111 

16 

870.24830 

6387.4618 

66.923305 

584176.57 

11.738844 

18* 

992.92178 

8071.9079 

64.713429 

752652.80 

18.598280 

22 

1244.3838 

12031.649 

63.844218 

1161^52.4 

39.663503 

30 

1705.8846 

2213.8303 

58.869145 

2167300.1 

180 . 90261 

38 

2151.1125 

34644.219 

54.543873 

3428287.3 

425.59583 

46 

2721.8130 

49616.932 

50.715587 

5300510.7 

647.48641 

54 

3476.8072 

67782.124 

47.445993 

8224919.4 

808.12281 

62 

4427.3496 

89914.139 

44.721666 

12693614.0 

924.52635 

70 

5617.1983 

116914.58 

42.479115 

19538075.0 

994.29298 

78 

7100.7933 

149892.68 

40 . 646909 

30033288.0 

1026.7656 

86 

8964.8019 

190264.90 

39.157532 

46305429.0 

1040.6460 

94 

11391.601 

240036.72 

37.954541 

72607238.0 

1045 . 6924 

102 

14785.568 

302500.64 

36.996442 

119039180 

1047.0206 

108 

18585.407 

361322.79 

36.152776 

184333900 

1047 . 0806 

*  End  of  tilt  phase, 
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TABLE  V   (Continued) 


c. 

n±   =   1.5 

Ty  -  16  sec 

U     =  18° 
m 

T 

V 

Y 

tf 

E 

VD 

sec 

fps 

ft 

deB 

ft-lb/slug 

fps 

k 

67.000877 

132.68853 

89.999999 

6513.6862 

8 

137 . 9621+2 

5I+I.2I+719 

89.999999 

26930.929 

.060251961 

12 

212.91698 

121+1.5799 

89.999999 

62613 .  kjk 

.3177309^ 

16 

291.88676 

22^9.6993 

89.999999 

111+980.88 

.9301+75 1+6 

20 

37^.09575 

3580.1953 

88.110768 

185163.20 

2 . 3827188 

2k 

1+59.5^522 

5236.3599 

82.72l)-083 

27^65.82 

1+.  8911670 

28 

551.81331 

7208.5201 

76.258076 

381+176.26 

8.2929821+ 

31* 

626.5981+2 

8900.2976 

72.77^338 

1+82671.1+1 

11.335807 

35 

732.62606 

lli+5l+.l+3^ 

71.827600 

636906.OO 

15 .  98271+5 

^3 

961k  122^ 

17603.6I+8 

65.55221^5 

103111+6.7 

31.150525 

51 

1207.21+51 

250I+5 . 312 

59.270032 

153^529.5 

72.761+797 

59 

11+50,8226 

33523.212 

53.090988 

2131020.6 

157.13915 

67 

1706.2258 

1^2720.706 

11-7.055876 

2830101.1+ 

278.08037 

75 

2027.7^50 

52525.718 

1+1. 291939 

37^5  81+0.0 

386.30959 

83 

21+23.2376 

62935 • ^90 

J5.98500I+ 

I+960929.9 

1+78. 6681+1 

91 

2898.1+151 

73921.632 

31. 228^59 

6578763.I+ 

55^.25915 

99 

31+52.0813 

851+22.589 

27.038602 

8706823.I 

619.8I+83I+ 

107 

1+089. 971+8 

973^9.100 

23.3800U9 

Hl+96062.0 

676.75^58 

115 

1+821.1+522 

109613.97 

20. 200?  1+0 

1511+9926.0 

72I+.  967^0 

123 

5658.59^ 

122137.53 

l-(,kkkOkl 

19939505.0 

765.06602 

131 

662I.I699U 

13^856.15 

15.057^69 

26262321.O 

79^.20610 

135 

715^.6835 

1I+1271.78 

13.987885 

3011+0033.0 

807.78328 

ll+7 

9009.8838 

160685.28 

11.208726 

I+5758899.O 

8I+I+.95I+17 

155 

10523.81^8 

173759.37 

9.6691722 

6096623 3.0 

867.32021 

163 

1235^.673 

186986.81 

8.3^6301+9 

82335093.0 

888.652I+7 

171 

1^652.533 

2001+90.62 

7.21791+99 

11379896.0 

910.07619 

180 

I8l8l.017 

216305. 1+7 

5.8239900 

17223I+II.O 

936.87I+76 

*  End  of  tilt  phase. 
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TABLE  V  (Continued) 


D. 

H--U 

I 

T     -  16-  sec 

12* 

u     =  ±8? 
m 

T 

V 

Y 

1 

E 

VD 

sec 

fps 

ft 

dep; 

ft-lb/slug 

fps 

k 

132.98571 

262.85677 

86.01*007 

17299.766 

.051083963 

8* 

Zjk.5kl39 

1061.2315 

78.231773 

71832.250 

.537^8828 

12 

1*25.70506 

2ko6.96kl 

76.811880 

16805 1*.  18 

1.551*81*19 

20 

755.87968 

678I*.  91*1*9 

69.I9187I* 

503976.21 

8.1*062369 

28 

1118,9996 

13^04.599 

62.519935 

1057360.3 

31.113581* 

36 

11*73.1912 

2201*8.882 

56.628937 

1791*51*8.0 

117.69071* 

IA 

1797.8701 

32228.993 

51.2W01 

2653105.7 

295.59529 

52 

2l81*.0l*12 

1*3659.1*80 

1*6.21*2080 

3789720.3 

1*81.81*897 

60 

2687.623^ 

56560.978 

1*1.736271* 

5^31^55.5 

629-6861*3 

68 

3312.1201* 

7120^.977 

37.809023 

7776023.5 

71*6.65796 

76 

i*06l*.6675 

87822.151 

31*.  1*1*1*157 

IIO86355.O 

81*0.0391*6 

81+ 

1*967.0819 

106668.91 

31.587081* 

15767923.0 

907.7^155 

92 

601*5.2513 

128059. k2 

29.17^501* 

22392722.0 

951.86607 

100 

7336.8l6l 

152389.86 

27.11*1*237 

318171*3^.0 

976.00713 

108 

889^.5880 

180171.00 

25.1*1*0375 

1^5353679.0 

991.09169 

116 

10821.625 

212127.15 

2k. OI6065 

65378775.0 

999.1*0505 

12k 

13307. ^39 

21*9378.99 

22.835027 

965671*97.0 

1003.5360 

132 

16759.873 

293825 • 31 

21.51*6871* 

11*9900220 . 

IOO5.0832 

135 

l8l*81*.668 

312997.36 

21. 25093  It 

I809II87O. 

1005 . 2250 

*  End  of  tilt  phase. 
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CHAPTER  7 


DISCUSSION  OF  RESULTS 


The  status  of  the  vehicle  at  burnout  is  of  prime  importance. 
This  vehicle  status  is  best  described  by  the  burnout  quantities,  V,  ,  Y^,  E^, 
and  &  -k,  as  displayed  in  Figs.  5  through  18. 

These  burnout  values  are  not  shown  for  the  vehicle  where  the  value 
of  n.  is  1.5  and  the  Ty  is  1  second,  due  to  the  fact  that  the  vehicle  passes 
through  the  horizontal  and  heads  back  towards  the  earth  before  reaching 
burnout.   This  result  is  readily  explained  since  the  vehicle  is  relatively 
slow  and  turns  at  a  low  altitude  where  drag  losses  are  very  large. 

The  vehicle  does  reach  burnout  conditions  for  the  other  values  of 

Ty,  however.   For  the  case  where  n-  is  1.5>  the  value  of  V,  peaks  at  a  very 

low  value  of  U  .   The  larger  values  of  U  cause  the  relatively  slow  vehicle 
m         &  m  J 

to  turn  more  at  a  low  altitude.   This  reduces  V,  due  to  the  fact  that  the 

b 

vehicle  operates  for  a  longer  period  in  dense  atmosphere  where  the  drag 

velocity  loss  is  very  large. 

As  T  is  increased,  the  maximum  value  of  V,  obtained  for  an  n.  of 
v  *  b  i 

1.5  occurs  at  the  higher  values  of  U  .   This  happens  as  a  result  of  the  in- 
creased velocity  and  the  higher  altitude  reached  before  the  tilting  is 
commenced.   Drag  velocity  loss  is  much  less  under  these  circumstances.   It 
is  true  that  gravity  velocity  loss  increases  as  T  is  increased  but  it  does 
not  offset  the  reduction  of  the  velocity  loss  due  to  drag.   This  fact  is 
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further  borne  out  by  observing  that  the  curves  indicate  higher  burnout  alti- 
tudes are  reached  as  the  value  of  U  is  reduced. 

m 

The  burnout  altitude  reached  varies  inversely  with  the  value  of  Um 
and  directly  with  the  value  of  T  for  the  vehicle  with  an  n.  of  1.5-   Both 
of  these  phenomena  are  readily  explained  since  the  altitude  attained  is  a 
direct  function  of  the  vertical  component  of  the  velocity  vector  which  is 
directly  affected  by  these  two  factors. 

The  value  of  tf  ,  of  the  vehicle  for  the  values  of  n.  »  1.5  generally 

decreases  with  increased  values  of  IT,  and  with  an  increase  in  T  .   A  large 

m  v 

U  allows  a  greater  amount  of  turn  of  the  vehicle,  hence  a  smaller  Hi  ■,.   At 

the  higher  values  of  T  the  vehicle  has  a  greater  velocity  before  commencing 

the  turn,  and  hence  does  not  get  turned  as  much  since  the  gravity  vector 

causing  the  turn  after  U  is  reached  is  small  in  comparison  to  the  vertical 

m 

component  of  the  velocity  vector  of  the  vehicle. 

Generally,  the  burnout  energy  of  the  vehicle  follows  the  trend  of 
the  velocity  at  burnout.   This  is  explained  by  the  fact  that  energy  is  directly 
proportional  to  the  square  of  the  velocity. 

As  the  value  of  n^  is  increased  the  burnout  velocity  attained  be- 
comes less  and  less  dependent  upon  the  value  of  Um,  which  is  indicated  by  the 

fact  that  the  curves  of  V-u  versus  U  tend  to  flatten  out  as  the  value  of  n. 

D        m  x 

increases.   Vertical  flight  time  does  not  affect  the  burnout  velocity  to  any 
great  degree  and  at  higher  vertical  flight  times  the  value  of  U  appears  to 
have  less  affect  on  the  V-^  reached.   The  basic  reason  for  these  effects  is 
that  with  a  high  n.  and  T  the  vehicle  is  out  of  the  very  dense  atmosphere 
before  turning,  and  hence  the  drag  velocity  loss  is  relatively  small. 

As  n.  is  increased  the  burnout  altitudes  become  less  and  less  at 

i 

the  low  values  of  Um  and  just  slightly  more  at  the  high  values  of  U .   At 
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the  higher  values  of  U  the  high  n.  vehicle  does  not  get  turned  as  much  as 
o  m       °   1 

the  lower  n.  vehicle  and  hence  attains  a  slightly  higher  burnout  altitude. 

Vertical  flight  time  seems  to  have  little  effect  on  the  burnout 
velocity  attained.   It  appears  to  be  a  slight  factor  at  low  values  of  T  , 
but  for  values  of  Ty  greater  than  8  seconds,  vertical  flight  time  has  no 
appreciable  effect  upon  the  burnout  velocity  reached.   The  value  of  "tf .  de- 
creases with  an  increase  in  n.  for  the  reasons  previously  explained. 

The  drag  velocity  loss  of  the  missile  for  the  various  trajectories 
studied  is  shown  in  Figs.  2k   through  27.   The  variation  of  this  loss  with  #  ^ 
and  n^  for  four  values  of  T  is  shown.  The  more  prominent  indications  of 
these  results  are: 

(1)  High  nj_  missiles  have  the  highest  drag  velocity  loss. 

(2)  The  drag  velocity  losses  decrease  slightly  with  an  increase 
in  vertical  flight  time. 

(3)  Drag  velocity  losses  greatly  increase  as  the  missile  attitude 
approaches  the  horizontal  at  burnout. 

The  higher  drag  velocity  losses  accompanying  an  increase  in  the 
value  of  n.  is  a  result  of  the  higher  velocities  attained  by  the  vehicle  at 
lower  altitudes.   Since  drag  force  is  directly  proportional  to  the  square  of 
the  velocity  and  to  the  density  of  the  atmosphere,  this  force,  and  hence  the 
resulting  velocity  loss,  is  large  for  a  high  n.  missile.   Since  the  drag 
velocity  loss  may  be  expressed  as: 

V»  =  So  /§  * 

it  must  be  realized  that  a  slow,  large  vehicle  would  experience  a  lower  drag 
velocity  loss  than  would  a  smaller,  faster  missile. 
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The  drag  velocity  loss  decreases  slightly  with  an  increase  in  ver- 
tical flight  time  because  the  missile  is  at  a  higher  altitude  and  hence  in 
a  less  dense  atmosphere  "before  it  commences  to  turn.   It  therefore  spends 
less  time  in  the  denser  atmosphere  of  low  altitudes.  This  decrease  in  drag 
velocity  loss  as  indicated  in  the  aforementioned  plots,  is  not  as  great  as 
would  be  expected.   The  drag  loss  is  plotted  versus  #  ,  ,  and  to  get  to  the 
same  value  of  ^-^   for  a  high  vertical  flight  time  as  for  a  low  vertical 
flight  time,  the  vehicle  must  be  turning  at  an  angle  of  attack  for  a  longer 
period  of  time  since  the  missile  has  a  greater  velocity  at  the  start  of  the 
tilting  phase.   This  increased  time  of  tilt  with  an  angle  of  attack  increases 
the  drag  coefficient  due  to  the  induced  drag  present  while  this  condition 
exists. 

Nearly  the  same  reasoning  applies  to  the  condition  of  increased 
drag  velocity  loss  for  a  smaller  value  of  ^  attained  for  a  particular 
vehicle  with  a  given  n.  at  a  certain  T  .   The  smaller  value  of  7$",  requires 
that  the  missile  be  tilted  at  an  angle  of  attack  for  a  longer  period  and 
hence  the  drag  coefficient  is  again  increased. 

The  values  of  n-r  which  are  plotted  against  U  in  Figs.  28  through 
31  show  the  maximum  negative  lift  load  factor  encountered  during  a  trajectory 
of  specified  Tv  and  Um.   Lift  is  always  negative  in  the  tilt  phase,  that  is, 
in  the  direction  of  rotation  of  the  vehicle.   Lift  load  factor  increases 
rapidly  with  increased  U  for  the  high  n.  trajectories.   The  slope  of  the 
lift  load  factor  curve  increases  as  T  is  increased,  as  could  be  expected 
from  the  higher  dynamic  pressures  caused  by  the  higher  velocities  associated 
with  long  vertical  flight  times.  An  interesting  aspect  of  the  higher  T 
plots  is  that  the  maximum  nT  seems  to  level  off,  reaching  a  maximum  of  about 
1.2  g's  for  an  n.  of  3.0.   In  these  cases  the  missile  has  reached  the  less 
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dense  portions  of  the  atmosphere,  where  the  extremely  low  density  has  offset 
the  increased  velocity,  and  the  maximum  n  encountered  has  become  essentially 

Li 

constant.  The  large  values  of  nL  encountered  for  high  n^  trajectories  and 
the  accompanying  bending  moments  are  too  great  for  most  contemporary  liquid 
fueled  vehicles.   Use  of  trajectories  of  this  nature  require  the  heavier 
structural  design  associated  with  solid  fueled  vehicles. 

It  might  be  mentioned  at  this  point  that  the  rate  of  tilt  should 
have  a  considerable  affect  on  lift.   This  program  uses  a  tilting  rate  of  two 
degrees  per  second,  which  is  considered  a  nominal  rate  for  a  large  rocket- 
powered  vehicle  control  system  to  achieve,  but  a  minimum  rate  to  reach  the 
maximum  programmed  tilt  angle  within  a  reasonable  time.   The  two  degree  per 
second  tilt  rate  was  selected  to  keep  the  lift  load  factor  within  practical 
bounds.   A  study  of  methods  for  obtaining  minimum  lift  loads  through  differ- 
ent types  of  tilting  programs  is  of  interest,  but  beyond  the  purposes  of 
this  paper. 

The  variation  of  tilt  time  necessary  to  reach  a  "tf,  of  30  for 
four  different  vertical  flight  times  is  plotted  against  n^  in  Fig.  27-   The 

curves  indicate  that  the  tilt  time  increases  as  both  T  and  n.  increases. 

v       1 

The  higher  the  value  of  T  used,  the  greater  the  velocity  of  the 
vehicle  becomes  before  the  vehicle  starts  to  turn.   Since  this  is  the  case, 
it  takes  much  longer  to  turn  the  vehicle  from  the  vertical  position  to  a  ^  , 
of  30  at  the  higher  value  of  T  than  at  the  lower  values.   This  same  type 
of  reasoning  may  also  be  applied  to  the  second  observation,  in  that  at  high 
values  of  n.  the  vehicle  gains  greater  velocities  sooner  than  at  low  values 
of  n . .   It  therefore  requires  a  greater  tilt  time  to  reach  a  ^ ,  of  30°  at 
the  higher  values  of  n. . 
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The  optimization  program  undertaken  in  this  paper  is  based  on  max- 
imum specific  energy.   The  combination  of  T  and  U  which  would  give  maximum 
burnout  energy  for  a  specified  burnout  angle  are  determined.   Energy  is  max- 
imized in  this  manner  for  the  three  higher  values  of  n. .   The  resulting  values 
of  E^,  U  and  Tv  are  plotted  versus  burnout  angle  in  Figs.  19  through  21. 
It  is  not  possible  to  obtain  a  good  optimization  for  an  n.  of  1.5  with  the 
data  available.  At  this  value  of  n-^,  burnout  energy  is  extremely  sensitive 
to  Ty  and  U ,  and  any  attempt  at  optimization  requires  a  large  number  of  tra- 
jectories.  Accordingly,  no  optimization  is  included  for  an  n.  of  1.5 • 

Good  maximum  energy  points  are  found  for  burnout  angles  below  about 

thirty  degrees.   In  the  range  of  thirty  to  fifty  degrees,  two  approximately 

equal  maximum  energy  points  appear.   One  of  these  points  occurs  at  a  value 

of  T  consistent  with  the  maximum  points  found  in  the  thirty  degree  and  below 

range,  while  the  other  occurs  at  the  minimum  T  ,  which  is  one  second.  Above 

the  thirty  to  fifty  degree  range  the  maximum  energy  points  are  found  at  the 

minimum  values  of  T  and  U  for  the  particular  burnout  angle.   These  same 

v      m        ^  ° 

characteristics,  in  varying  degree,  are  found  for  each  value  of  n. .   The 
overlapping  in  the  plots  of  Tv  and  Um  versus  burnout  angle  in  the  thirty  to 
fifty  degree  range  show  that  in  this  area  maximum  energy  can  be  obtained  by 
using  either  of  two  combinations  of  Tv  and  U  .  This  effect  is  undoubtedly 
caused  by  the  non- linear  action  of  drag  on  the  trajectories.  A  point  is  reached 
for  each  n.  where  the  beneficial  effects  of  early  tilting,  and  consequent  earl- 
ier alignment  of  thrust  and  velocity  vectors  in  the  desired  direction,  are 
overcome  by  the  higher  drag  losses  associated  with  large  programmed  tilt  angles 
at  low  altitudes.  At  this  point  it  is  necessary  to  use  a  period  of  vertical 
flight  time  to  get  the  vehicle  out  of  the  denser  portions  of  the  atmosphere 
during  the  tilting  maneuver. 
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The  amount  of  vertical  flight  time  required  to  maximize  burnout 
energy  for  low  "burnout  angles  does  not  increase  indefinitely  as  burnout  angle 
approaches  zero,  but  tends  to  level  off  in  the  neighborhood  of  l6  to  20  seconds. 
This  is  especially  evident  in  Fig.  21  where  n*   is  3.0.  The  high  velocity 
reached  at  the  end  of  a  long  vertical  flight  time  causes  an  increased  negative 
lift  load  during  the  tilting  phase,  which  increases  the  drag  coefficient, 
thereby  reducing  burnout  velocity  and  burnout  energy. 

It  is  interesting  to  note  that  at  burnout  angles  of  about  ten  degrees 
and  below,  the  maximum  burnout  energies  for  all  three  values  of  n.  fall  very 
close  together.   This  indicates  that  for  very  low  burnout  angles  the  advantage 
of  a  high  initial  acceleration  rocket  is  questionable,  since  the  same  burnout 
energy  can  be  obtained  using  a  lower  initial  acceleration  with  lower  aerodynamic 
loads.   At  other  values  of  burnout  angle  the  higher  burnout  energies  obtained 
from  high  n.  rockets  are  apparent.  The  point  at  which  maximum  overall  burnout 
energy  is  reached  starts  at  a  burnout  angle  in  the  vicinity  of  twenty-five 
degrees  for  an  n.  of  2.0,  and  moves  in  the  direction  of  increasing  burnout 
angle  as  n.  increases.   In  this  case,  the  increased  drag  associated  with  high 
initial  acceleration  and  low  burnout  angle  causes  the  maximum  overall  energy 
points  to  fall  at  higher  burnout  angles  for  the  higher  values  of  n. . 

Values  of  burnout  velocity  and  altitude  for  the  maximum  energy  condi- 
tions discussed  before  are  shown  in  Fig.  22.   It  logically  follows  that,  since 
energy  is  a  function  of  velocity  and  altitude,  these  curves  are  coincident  in 
the  same  manner  as  the  maximum  energy  curves  at  low  burnout  angles.  At  higher 
burnout  angles,  burnout  velocity  increases  and  burnout  altitude  decreases  as 
n.  is  increased. 

Each  representative  trajectory  tabulated  in  Table  V  is  the  nearest 
trajectory  available  to  the  overall  maximum  energy  case  for  that  particular  n . . 
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This  is  the  optimum  trajectory  for  the  burnout  angle  listed  only,  and  it  can- 
not be  said  that  it  is  the  optimum  trajectory  for  any  other  attitude  reached 
before  burnout.   Naturally  it  is  not  applicable  to  a  burnout  angle  lower  than 
that  listed.   The  question  arises  whether  or  not  more  energy  is  obtained  by 
using  the  tabulated  trajectory  for  the  overall  maximum  energy  case  until  the 
desired  attitude  is  obtained,  followed  by  constant  attitude  thrust  until 
burnout,  than  by  using  the  maximum  energy  trajectory  for  the  burnout  angle 
corresponding  to  the  attitude  desired.   It  seems  that  if  the  vehicle  is  above 
the  denser  atmosphere,  the  energy  generated  after  constant  attitude  thrust 
is  started  would  be  about  the  same  as  that  for  the  gravity  turn.   A  constant 
attitude  thrust  program  would  give  higher  altitude  with  lower  velocity  than 
the  gravity  turn,  although  the  actual  difference  between  the  two  programs 
would  depend  upon  the  time  of  application  of  the  constant  attitude  thrust 
program.   If  the  vehicle  reaches  the  desired  attitude  in  the  early  tilt  phase, 
before  leaving  the  denser  atmosphere,  a  constant  attitude  thrust  program 
would  involve  lift  loads  and  additional  drag,  but  it  would  get  the  vehicle 
out  of  the  sensible  atmosphere  sooner.   In  either  case  the  non-linearity  of 
the  problem  would  require  a  separate  computation  for  each  situation. 
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CHAPTER  8 


CONCLUSIONS 


Determining  the  optimum  powered  flight  trajectory  for  a  large 
rocket  powered  vehicle  is  a  complex  problem,  which  is  strongly  influenced 
by  the  desired  trajectory  burnout  angle  and  the  rocket  initial  thrust-to- 
weight  ratio.   The  burnout  angle  may  be  considered  a  design  parameter  for 
the  booster  trajectory,  since  different  burnout  angles  are  required  for  dif- 
ferent missions.   This  paper  shows  that  a  combination  of  vertical  flight 
time  and  initial  tilt  angle  to  give  maximum  energy  at  burnout  can  be  deter- 
mined for  any  desired  burnout  angle.  There  is,  in  addition,  one  value  of 
burnout  angle,  which  gives  maximum  burnout  energy,  for  each  value  of  initial 
acceleration.   In  this  manner  an  optimum  booster  flight  trajectory  is  avail- 
able for  any'desired  burnout  angle. 

Usually  the  vehicle  with  the  higher  initial  acceleration  will  have 
the  greater  burnout  energy.   At  low  burnout  angles,  however,  in  the  zero  to 
ten  degree  range,  values  of  burnout  energy  for  a  wide  range  of  initial  accel- 
erations closely  coincide.   For  low  burnout  angles,  therefore,  the  initial 
acceleration  of  the  vehicle  is  of  little  consequence  with  respect  to  maximum 
energy  optimization.   In  fact,  it  can  be  said  that  a  lower  initial  accelera- 
tion is  preferable  for  low  burnout  angles,  since  the  high  lift  load  factors 
associated  with  high  initial  accelerations  are  avoided. 
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The  time  required  to  tilt  the  vehicle  from  the  vertical  to  a  point 
where  the  angle-of-attack  is  zero  is  also  quite  high  for  high  acceleration 
vehicles.  The  procedure  involving  a  relatively  small  initial  tilt  angle 
followed  by  a  gravity  turn  to  the  desired  burnout  angle  is  no  longer  feasible 
with  high  initial  accelerations,  which  leads  to  the  conclusion  that  the  tilt 
phase  for  high  initial  acceleration  vehicles  must  be  programmed  in  its 
entirety. 
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APPENDIX  A 


ATMOSPHERIC  DATA 


The  atmospheric  data  used  is  based  on  the  ARDC  model  atmosphere 
of  1959  as  described  in  Ref .  3«   In  order  to  facilitate  computer  procedures 
the  density  ratio  (  f/fi,)  and  sonic  speed  data  are  treated  in  a  simplified 
manner . 

Utilizing  the  atmospheric  density  data  of  Ref.  3  a  plot  of  the 
ratio  of  £  /q     is  made  extending  from  sea  level  to  an  altitude  of  ^00,000 
feet,  as  shown  in  Fig.  32.  The  resulting  curve  is  divided  into  four  seg- 
ments which  are  accurately  approximated  by  appropriate  exponential  functions, 
The  resulting  segments  of  the  curve  and  respective  describing  exponential 
functions  representing  the  density  ratios  are  shown  in  Table  VI. 

Sonic  speed  is  plotted  versus  altitude  from  sea  level  to  an  alti- 
tude of  U00,000  feet  as  shown  in  Fig.  33-   The  curve  results  in  a  series  of 
five  straight  line  segments.   Straight  line  functions  are  used  to  describe 
these  segments  of  the  curve.  Table  VI  displays  these  functions  and  their 
respective  areas  of  applicability. 
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TABLE  VI 


ATMOSPHERIC  DATA  APPROXIMATION  FORMULAS 


1 

Altitude 
(ft) 

Density  Ratio 

, 

Speed  of  Sound 
(fps) 

0 
36,800 

e-Y/32,000 

1120-.001H7  Y 

82,500 

1.65  e-Y/20,800 

968.08 

120,000 

1.65  e-Y/20,800 

813.78      .00]87  Y 

168,000 

§51  e-Y/25,200 

813.78      .00187  Y 

263,000 

.51  e"V25,200 

1625 -.00298  Y 

i+50,000 

?T  e-Y/l7,000 

81*6.5 

7^ 


800 


900        1000        1100 
Speed  of  Sound  (fps) 


1200 


Fig.  33   Variation  of  speed  of  sound  with  altitude. 
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